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We consider a variant of the entanglement of assistance, as independently introduced 
by D.P. DiVincenzo et al. (quant -ph/9803033) and O. Cohen (Phys. Rev. Lett. 80, 
2493 (1998)). Instead of considering three-party states in which one of the parties, the 
assistant, performs a measurement such that the remaining two parties are left with 
on average as much entanglement as possible, we consider four-party states where two 
parties play the role of assistants. We answer several questions that arise naturally 
in this scenario, such as (i) how much more entanglement can be produced when the 
assistants are allowed to perform joint measurements, (ii) for what type of states are 
local measurements sufficient, (iii) is it necessary for the second assistant to know the 
measurement outcome of the first, and (iv) are projective measurements sufficient or arc 
more general POVMs needed? 
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1. Introduction 

Consider a pure state of four qubits located at parties A, B, C, and D. Whether this state is 
entangled or not, the reduced density matrix of the AB system, /?ab, might only possess very 
little entanglement, or it might even be separable. But the parties C and D, both knowing the 
full state, can help increasing this entanglement by performing measurements on their qubits 
and communicating their results to A and B. This leads to the concept of the entanglement 
of assistance [1] , defined as the maximum possible average entanglement the assistant parties 
A and B can create between C and D. This was originally introduced by Cohen [2] in the 
notion of maximal hidden entanglement and later independently rediscovered by the authors 
of Rcf. [1]. It follows from a result of Hughston, Jozsa, and Wootters (HJW) [3] that if 
the assistant parties are allowed to perform joint measurements, then the entanglement of 
assistance depends only on pab, and as such it is considered in Ref. [1] as a measure of 
two-party entanglement. We denote this quantity by Here we ask the question, what is 

a We use a musical notation where t\ ('natural') signifies that only local measurements (and classical communi- 
cation) are included, whereas fl ('sharp') signifies that joint measurements are included. Later on we will also 
use b ('flat') to signify that only local projective measurements are included. 
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the entanglement of assistance if the assistant parties are restricted to performing only local 
measurements? This quantity, which we denote E\ is not intended as a new measure of 
two-party entanglement, indeed it will in general depend on the full four-party state. Clearly 

< E\ but can this upper bound be reached, and if so for what states? 

Before attacking the problem we look at some simple examples to catch a glimpse of the 
properties of E* in relation to EK For example consider the four-party GHZ state, 

_ 10000) + 

IWabcd ^ • (1) 

If C and D both measure in the basis {|+), |— )}, defined as |±) = (|0) ± |1)), A and B will 
always end up with a maximally entangled state, so that in this case the two quantities are 
equal, E* = E^ = l b On the other hand, for the state c 

|V-)abcd = |$ + )ac|$ + )bd, (2) 

no local measurements by C and D can change the fact that A and B are separable, and hence 
E^ = 0. But with a joint measurement in the Bell state basis, C and D can ensure that A 
and B end up in a maximally entangled state, as illustrated in Fig. 1, so that E$ = 1 for this 
state. This is obviously the largest difference we can get between E^ and E^. 




Fig. 1. A simple example of a state where joint measurements can create entanglement across 
a separable cut where local measurements cannot, is the entanglement swapping indicated here. 
C and D make a joint Bell state measurement on the initial state, Eq. (2), thereby creating a 
maximally entangled state between A and B. For this state the entanglement of local assistance 
is zero. 

A result which is useful for a certain class of states is the one by Walgate et al. [4] , stating 
that two orthogonal states of two qubits can always be locally distinguished. This result 
can be applied to immediately determine a set of states for which local measurements are 
sufficient. Indeed for the states where the optimal joint measurement is a projection onto at 
most two orthogonal states (with nonzero probability) we have E^ = EK A particular local 
measurement that achieves this can be found in Rcf . [4] . For this measurement we note that 
classical communication of the outcome of C's measurement to D is explicitly needed. 

Yet another example of a state for which communication between the assistants seems to 
be necessary is 

|^)abcd = J(|3> + )ab|00) C d + |3»~)ab|01) C d + |* + )ab|1+)cd + |*~)ab|1->cd) • (3) 

6 Here we use the standard measure of entanglement, E(\ip)) = — Trp^ log 2 Pa> see Ref. [10]. 

c Here and throughout we have denoted the four maximally entangled Bell states \^^) = -j= (|00) ± 1 1 1) ) and 

l* ± > = ^(|01)±|10)). 
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For this state one optimal local measurement is clearly the one where C measures in the 
{|0), | l)}-basis and D measures either in the {|0), |l)}-basis or in the {|+), |— )}-basis de- 
pending on the outcome of C's measurement. However, we will show in Sec. that, rather 
surprisingly, classical communication of measurement outcomes between the assistant parties 
is never necessary, as long as they restrict themselves to von Neumann measurements. In 
particular we note that for the above state, Eq. (3), C can measure in the {|0), |l)}-basis and 
D in the {-^+|[I1 5 Mz|il> j.^asis irrespective of C's measurement outcome. 

We note that it might be necessary for C and D to communicate to decide on which 
optimal basis to use, if there is no unique such basis. Moreover classical communication from 
CD to AB will of course always be necessary. Otherwise, from their point of view, AB would 
be stuck in the state pab no matter what is happening at C and D. 

For convenience we will take the concurrence of Wootters [5, 6] as our measure of entan- 
glement. Concurrence is a well-established measure of entanglement, in particular it is an 
entanglement monotone [7]. Moreover, the concurrence features in some elegant inequalities 
[8] , see also [9] , that are not satisfied by the entanglement of formation. 

Thus we introduce the concurrence of assistance. When we allow joint measurements we 
will denote this quantity by CK When we allow only local measurements we will distinguish 
two quantities, namely C b where we allow only local von Neumann measurements, and & 
where we allow more general POVMs. In the main part of the present paper we will restrict 
the assistant parties to local von Neumann measurements and classical communication. This 
as well as a scrutiny of how much freedom we have in choosing optimal (local or joint) 
measurements for a given state, will be the topic of Sec. . In Sec. we will show that classical 
communication of measurement outcomes between the assistants is not necessary. In Sec. we 
will find a decomposition that allows us to quickly determine (i) the concurrence of assistance, 
(ii) whether there is a local measurement with which we can achieve this, as well as (iii) the 
explicit form of this local measurement. In Sec. we will present some numerical results 
showing the advantage of joint measurements on arbitrary states, and in Sec. we discuss the 
effect of non-pro jective measurements and some possible extensions of the present work. 

2. Concurrence of Assistance 

The entanglement of assistance has been considered by DiVincenzo et al., Ref. [1]. Starting 
from a three party pure state \ip) abc> it quantifies how much pure state entanglement party 
C can make between A and B by doing local measurements on his part of the state. 

The reduced state pab can be described by infinitely many different pure state ensembles, 
£ = {pi, \(f>i)} all satisfying /?ab = X^i Pi I </>*)(</>» I > an d from each ensemble we can construct a 
purification as 



where {|«)c} is an orthonormal basis of C's system, which might have been extended by an 
ancilla system to increase its dimension. Only one of these purifications is the state |^>)abCi 
but the HJW-theorem [3], mentioned above, states that we can transform purifications cor- 
responding to different ensembles into each other by a local unitary operation acting on the 
system of the assistant party C alone, which means that, in particular, all the states \tps) of 
Eq. (4) can be reached from our state |^)abc by local operations on this system. 




(4) 
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A local measurement by C in the basis {|i)c} will make the AB-system end up in the 
state \cj)i) with probability pi, thus realizing £ as a classical statistical ensemble. Conversely 
any measurement that C can do on the state, will correspond to such an ensemble. These 
considerations led the authors of Ref. [1] to define the entanglement of assistance as 

£» = rnax^£(|<^)) . (5) 

i 

The measure of entanglement, E(\4>i)), used here is the von Neumann entropy of the marginal 
density matrix px — TrB|</>j)(0j|, coinciding for pure states with such other measures as the 
entanglement of formation and the entanglement of distillation, see e.g. Ref. [10]. 

We see that the entanglement of assistance is a function of the reduced state pab only, 
and as such it is suggested in [1] as a measure of two-party entanglement. An explicit formula 
has not been found, but for the case where the systems of A and B arc qubits, it was shown 
in [1] that < F(p, p), where p — pab, and where the fidelity function is defined as 

F(p,a) = TrvV /2 er/o 1/2 . (6) 

Here the Wootters tilde, Refs. [6, 5], is defined for pure states and mixed states respectively 
as 

W) = i°y®Oy)W), 

P = {<Jy <g> Gy) P* (<Ty <g> Gy) , (7) 

with the complex conjugation taken in the standard basis. 

In these two references, see also Ref. [11], an alternative measure of entanglement for two 
qubits, namely the concurrence, is defined. For pure states \tp) = a\00) + b\01) + c|10) + d\ll) 
this is 

C(\TP)) = Mi,)\=2\ad-bc\. (8) 

The main advantage of this measure is its simple algebraic form, which enables an analytic 
treatment of some of the problems to be presented in this paper. In analogy to Eq. (5) we 
therefore define the concurrence of assistance 4 as 

C»= max . (9) 

i 

The HJW-theorem ensures also in this case dependence on pab only. 

One specific ensemble describing our two-qubit state p = pab is the one found by Wootters 
in Ref. [6], £yy — {Qi, \ x i)} satisfying ^J<Mj{ x i\ x j) — where A, denotes the eigenvalues of 

the matrix \J p 1 / 2 pp 1 / 2 . A measurement yielding these |a;»)'s with probabilities will result 
in an average concurrence of 

C{£w) = ^2pi\{xi\xi)\ = ^2 K = F(pab,Pab) ■ (10) 

i i 

d We remark that the function E(C) relating the concurrence to the entanglement measure used in Eq. (5), is 
not linear [6], so in particular _E" ^ _B(C"), and hence the concurrence of assistance is really an alternative 
measure to the entanglement of assistance, rather than a simple function thereof. 
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The strong concavity property of the fidelity function, Theorem 9.7 of Ref. [12], tells us that 
for any ensemble £ = {71^, \4>i}} describing pab, this is an upper bound, 

c{£) = 5>iK&l&>l = I> F (l&X&l- l&X&l) 

^F^Tril^X^I.^TTil^X^I^ =F(pab,Pab), (11) 

so the value (10) reached with the ensemble £w is in fact equal to the concurrence of assistance, 

C s = F( Pab ,pab), (12) 

as was pointed out by C.A. Fuchs [13]. 

In the following we will focus on pure states of four qubits. For such states a new problem 
arises, namely the one of local assistance. How much concurrence can the assistant parties 
C and D make for A and B after local measurements on their own qubits, C\ and what 
measurement is the optimal one? These and other related questions will be the topic of the 
following sections. 

First we derive some general properties of the concurrence of assistance. From now on let 

|V) = ci|0000) + c 2 |0001> + • • • + ci 6 |llll) (13) 

be the initially shared pure state of four qubits located at parties A, B, C, and D respectively, 
and construct the matrix 

/ ci c 2 c 3 c 4 \ 

j£ _ C5 C6 Cj C 8 , . 

Cg ClO Cn Cl2 
\ C13 C14 C15 Ci 6 / 

By explicitly checking we see that 

Pab = AAt , (15) 

and by symmetry it then follows that pcd = X T (X T Y . A useful way to think of X is that if 
we write \tp) — |0oo)|OO) + • ■ • + |0n) 1 1 1) , then X is the matrix with the unnormalized |0»j)'s 
as column vectors, 

( \ I I I \ 

X = |0oo> |0oi> l^io) \4>n) ■ (16) 

VI I I I / 

A measurement done by C and D in their {|0), |l)}-bases with outcome |i)c |j)d will project 
the state of the AB-system onto \4>ij), and this outcome will occur with probability (4>ij\4>ij). 
The concurrence of this state, Eq. (8), is \(cj>ij |</>ij)|/(</>ij|</>ij), so the average concurrence 
obtained by this particular measurement is 

c = £l<4l<^)l- (i7) 

ij 

In terms of the matrix X introduced above this reduces to 

c = E|((K ^) x *) tx ) fc J = EK xT K ^) x ) fe J • ( 18 ) 
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Transformations of the type 

X — >XV, (19) 

where V is a (rectangular) right unitary matrix, VV^ = 1, do not change pab (15), and this 
is indeed the only type of transformation of X that leaves pab unchanged. This tells us the 
freedom we have in choosing an X satisfying Eq. (15) for a given state, which is completely 
equivalent to the freedom mentioned above (4) we have in choosing the ensemble. 

Suppose that V is square, and hence unitary (this corresponds to a von Neumann mea- 
surement). From Eq. (16) we see that the column vectors of XV are the final states of the 
AB system after C and D have done a (von Neumann) measurement in the standard basis 
on the state (XV) T {XV) T ^ = V T pcdV t \ or equivalently, a measurement on pcD in the 
basis consisting of the column vectors of V*. In general this involves joint actions on the CD 
system. 

The average concurrence of the AB-state after this measurement will be as in Eq. (18), 
but with X — ► XV, so the concurrence of assistance is the maximum 

C S = max]T|(y T QVU > ( 20 ) 

k 

where the matrix Q is defined as 

Q = X T {<7y ® <7y) X . (21) 

To find this maximum we apply some mathematical results, the first one being the ex- 
istence of Takagi's factorization, see Rcf. [14]. It states that if Q is a complex symmetric 
matrix we can find a unitary matrix U such that 

\ 

U r QU= " n " n a ° J , (22) 
er 4 / 

where tJfc are the singular values of the matrix Q. 

Our matrix Q of Eq. (21) is indeed symmetric, Qij — X)fcz(^ T )ifc( fT y ® o~ y )kiXij = 
J2ki(X T )ji(o~ y (g> o-y)ikXki — Qji, and together with the next mathematical result, namely 
that 

£|(V^Qy) fefc |<]>>, (23) 

k k 

for all right unitary matrices V, see App. 1, we get the following nice formula for the concur- 
rence of assistance, 

C» = £ \(U T QU) kk \ =^>. (24) 

k k 

Moreover we get a particular candidate for the optimal measurement, namely the von Neu- 
mann measurement projecting onto the column vectors of U* , where U is the unitary matrix 
found from the Takagi decomposition of Q, Eq. (22). 

This is, however, not the only optimal measurement. For our present problem it is espe- 
cially interesting to see whether there is an optimal measurement that is local. To see this 
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we need to specify what freedom we have in choosing matrices V that reach the maximum in 
Eq. (20). For convenience we shall assume that V is unitary, and hence restrict C and D to 
von Neumann measurements. 

Singular values are per definition real and positive, so from equation (24), C" = Ti(U T QU) 
= Tr(0 T U T QUO), showing that the unitary matrix UO, where O is an arbitrary real, ortho- 
gonal matrix (O t O = 1), is another choice of optimal measurement basis. Supposing that all 
the singular values of Q are nonzero, further multiplication by unitary matrices that does not 
change the sum of the absolute values in Eq. (24) , can only be with a diagonal unitary matrix, 
since these change only the phases of the entries. So in this case all unitary V's corresponding 
to an optimal measurement basis are of the form 



V = UO 



( 













\ 
































V 













) 



(25) 



where U is a particular matrix satisfying Eq. (22), O is a real, orthogonal matrix and 
ipi, . . . , ifi are arbitrary real numbers. 

3. Classical Communication is Not Necessary 

The results of the preceding section do not concern distinguishing local and joint measure- 
ments, but we are now in a position to treat this problem. For local measurements we have 
the possibility of communicating measurement outcomes between the assistant parties, and 
letting the outcome of one party's measurement determine the measurement basis to be used 
by the other party. As discussed in the introduction, one might expect this to be necessary for 
the optimal local measurement, but in this section we shall show that as long as the measure- 
ments are von Neumann measurements, classical communication of measurement outcomes 
between the parties is actually not necessary, no matter what state we are looking at. 

To show this, suppose that party C is the first party to do a measurement, and that he 
docs it in the basis W (by "a measurement in the basis W" we mean "a measurement in the 
basis consisting of the column vectors of W* v , as discussed above). Party D now has the 
opportunity to let her measurement basis depend on the outcome of C's measurement, so the 
unitary matrix V corresponding to their combined measurement will be of the form 

V = (W®1)( V ( ; °), (26) 

where W , V\, and V2 are all two-dimensional unitary matrices. We will now show, that among 
the best measurements that D can do following this measurement by C, there is one satisfying 
V\ = V2 , so that in total V = W ® V\ , and therefore D doesn't have to wait for the outcome 
of C's measurement. 

The concurrence of assistance of a state, with our restriction to local von Neumann mea- 
surements and classical communication, is according to Eq. (20) 

C» = m r -£\(V T QV) kk \ , (27) 
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where the maximum is taken over unitaries, V, either of the form (26) or of the similar form 
with the roles of C and D reversed, and where Q is defined by Eq. (21). Let 



(W T ® 1) Q {W ® 1) = ( Q 



* 

Q 2 



(28) 



where the *'s indicate that the off-diagonal blocks are irrelevant. After C's measurement the 
state has collapsed, so that either the first or the last two columns of the matrix X(W (g> 1) 
will be all zeros, see Eq. (16), and consequently the matrix in Eq. (28) will have become one 
of the two, 



Qi 




or 




Q 2 



(29) 



Since (W T <g> l) Q (W <g> 1) is complex symmetric, so are Qi and Q 2 , Eq. (28), and we can 
therefore write the Takagi-decomposition, Eq. (22), of these, Q\ — V*T,\V^ and Q 2 = 
V 2 *^2V 2 . From Eqs. (26), (27), and (28) the highest achievable average concurrence fol- 
lowing this measurement by C is 



m; 

Vi,V 2 



ax > 



V 2 T Q 2 V 2 



which, by the same argument as in Eq. (24), reduces to 

Tr(S!) + Tr(S 2 ) . 



(30) 



(31) 



Hence we see that V± and V 2 are one maximizing choice for V\ and V 2j and as in Eq. (25) we 
see that we will get the same average concurrence by choosing V\ and V 2 as 



Vi = V1O1 











and V 2 = V 2 2 











(32) 



where 0\, 2 are real orthogonal matrices and ct\, a 2 , (3i, and (3 2 are arbitrary real numbers. 
We want to show that we can find 0\ and 2 and a set of phases such that V\ = V 2 , or 
equivalently 



V^V 2 = Oi 











(33) 



But this is in fact always true, due to the following theorem: 

Given a unitary n x n matrix U we can find real orthogonal matrices 0\, 2 , so that 



U = Oi 



( e lSl 


V 







\ 









,i6n 



(34) 



/ 



The phases, which we will refer to as the phases ofU, Si (— -| < Si < J) are uniquely defined 
(up to permutations). The proof of this theorem, as well as how much freedom we have in 
choosing 0\ and 2l is given in Appendix B. 
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Hence it is always possible to fulfil Eq. (33), and therefore for C and D to make an 
optimal local measurement in a basis of the form V = W ® V\ . This means that classical 
communication of measurement outcomes is not necessary to improve the average concurrence 
achieved in this type of measurements, and in particular we note that it doesn't matter which 
one of the assistant parties, C and D, performs the measurement first. 

Having obtained this result we now also get a simple expression for C b . From Eqs. (27) - 
(31) it immediately follows that 

C b = max(Tr(S 1 ) + Tr(£ 2 )) , (35) 
w 

where W describes the measurement done by C, and Si and T, 2 are the singular value matrices 
of the two parts of Q after C's measurement, defined in Eq. (28). 

4. For What States are Local Measurements Sufficient? 

In this section we are going to find the set of states for which the assistant parties, C and D, 
have no advantage in using joint measurements as compared to local von Neumann measure- 
ments, i.e. C — CK Of course this automatically implies that also = CK 

4.1. rank(S) = 4 

Our starting point will be the expression, Eq. (20), for the concurrence of assistance. At 
first we assume that the initial state \ip) has rank(pAB) = 4, which is equivalent to assuming 
that all the singular values of the matrix Q of Eq. (21) are nonzero, i.e. rank(S) = 4. The 
remaining cases will be discussed separately below. For the states we consider here, the degree 
of freedom we have in choosing the optimal measurement basis V is completely determined by 
equation (25). Together with the result of the previous section that classical communication 
of measurement outcomes is not necessary for this type of measurements, this equation tells 
us that local von Neumann measurements can do as well as the best joint measurement if and 
only if we can find matrices Cq £ SO(4), D £ SU(4) (diagonal), and Ui,U 2 £ SU(2) such 
that 

UOiD = U 1 ®U 2 , (36) 

where U £ SU(4) is determined by Eq. (24). For convenience we are requiring all matrices to 
have determinant onef 

Now, define the unitary matrix T £ SU(4) by 



(37) 



This matrix has the useful property that for any two matrices U\,U2 £ SU(2) the matrix 
T (Ui ® U 2 ) T f £ 5*0(4), and conversely, for any matrix O £ SO (4) there exist U X ,U 2 £ SU{2) 
such that T^OT = Ux® U 2 , see Rcf. [15]. 

e This is a legal requirement, since it only amounts to multiplying the matrix with a constant phase factor. 
This does not change the concurrence, Eq. (20), and only changes the measurement basis (columns of V) by 
an overall phase factor. 
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Applying this result to the criterion Eq. (36), we see that it becomes 

TU = 2 TD^O\ , (38) 

where 2 = T {U\ <g> U2) . The matrix TU is unitary so we can apply the decomposition 
(34) to find its phases, i.e. to find 3 ,0 4 £ SO (4) and E £ St/ (4) (diagonal) such that 

TU = 3 EOj . (39) 

The criterion (38) now becomes: when can we find O5, 0% £ SO(4) and D £ SU (4) (diagonal) 
such that 

E = O^TD^Ol , (40) 

where E is determined from the state \tp) by Eq. (39). The matrix TD^ is also unitary, and 
accordingly the phases of that matrix can also be determined by writing it in the form of 
Eq. (34). When varying D there must be some constraints on what the phases of TD^ can be. 
(If all choices were possible, it would mean that local measurements were always sufficient, in 
contradiction with several counterexamples, as discussed in the introduction). 

By explicit calculation, following the lines of App. B, we find that with the diagonal entries 
of D equal to {e^ 1 , e iV2 , e iV3 , e^ 4 }, the phases of TD^ turn out to be 

[e J *,e- J *,e-<*-f),e<*-t)} , (41) 

where $ = (j(<Pi + <Pa — <f2 — V3) mod -|) lies in the interval < $ < -|. Hence our criterion 
for local von Neumann measurements to be sufficient is that there has to exist a number $ 
such that the phases of TU, Eq. (39), arc of the form (41). 
U is found from the state X according to (21), (22), 

U T X T (<Ty ® <7y) XU = £, (42) 

where S is a diagonal matrix containing the singular values of Q = X T (a y (8) a y ) X, which 
according to Eq. (24) gives us the concurrence of assistance, 

C tt = Tr (£) . (43) 

The matrix 

- 1 



2 



( l+i \ 

l + i 1 — i 

1 — i l + i 

\ -l + i l + i ) 



(44) 



is symmetric, and so is so Eq. (42) is of the form A T A = B T B, where A and B are 
nonsingular. This determines uniquely a complex orthogonal matrix 0, (CI = (BA~ 1 ) T ), 
satisfying A = QB, or 

^/Oy ® OyXU = ftVs . (45) 

Putting together this result and the decomposition of TU of Eq. (39), we end up with the 
following useful decomposition of the state, 

X = ^/a y <g> a y ^nVY,PiFP?T . (46) 
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Here we have defined the diagonal unitary matrix F = OjE^O?, and the orthogonal matrices 
Pi = O4O7 and P 2 = O3O7. O7 is an orthogonal matrix permuting the phases of E to ensure 
that P2 (and Pi) have determinant one, and that, for the states where the phases of F are as 
in Eq. (41), the ordering is either as in that equation or with the first and the fourth phase 
interchanged (only one of these orderings will be consistent with det P2 = det Pi = 1). 

In this decomposition S gives the concurrence of assistance as in Eq. (43), and F tells us 
whether it is possible to obtain this by local von Neumann measurements, according to the 
criterion we found, which was that its diagonal entries have to be the phase factors of Eq. (41) 
for some $. 

The product of the last four terms in this decomposition is W = PxFP^T ', which gives 
us information about which basis we should use in the optimal local measurement. Applying 
the mathematical result stated after Eq. (37) to the orthogonal matrix P2 = TUP\F 1 we get 
that there exist ?7i, f/2 € SU{2) such that 

P t P 2 T = Ui®U 2 . (47) 

But this is also equal to T^ P 2 T = UPiFT, and for the states where local measurements are 
indeed found to be sufficient, the diagonal entries of F are now either as in Eq. (41) or as in 
this equation but with the first and the fourth phase factor interchanged. In both cases we 
check explicitly that the matrix FT is of the form of a real orthogonal matrix times a diagonal 
unitary matrix. Hence, according to Eq. (25), the matrix P 2 T = UP\FT of Eq. (47) gives 
a measurement basis with which we can achieve the concurrence of assistance, and, most 
importantly, this measurement basis corresponds to a local measurement. 

Thus, in summary, writing the state of four qubits in the decomposition (46) provides a fast 
way to find its concurrence of assistance (£), to decide whether local projection measurements 
by C and D are sufficient to achieve this (F), and what those local measurements would be 
OP2). 

It is interesting to count the (continuous) degrees of freedom we have in choosing states 
for which local von Neumann measurements are sufficient, to see how much space these states 
take up in the full state space. Normalization of the state X implies that Tr (XX*) = 
Tr (Sl£f^) = 1, so the complex orthogonal matrix O has 12 — 1 = 11 degrees of freedom. The 
singular value matrix, s/Ti, has 4, the two real orthogonal matrices, Pi and P2 have each 6. 
The unimportance of the overall phase led us to assume det(P) = 1, so that in general F has 
3 degrees of freedom, thus counting up to the total number of 30 = 32 — 2 degrees of freedom 
for an arbitrary pure state of four qubits \ip) € C 16 . For states where local von Neumann 
measurements are sufficient, F has only one degree of freedom, namely the $ of Eq. (41), 
giving a total of 28 degrees of freedom in choosing these states. Thus, although these states 
span a set of measure zero on the full state space, they have only two continuous degrees of 
freedom less than the full state space. 

4.2. The rank (S) < 4-cases 

The above derivation was completely general, except for one point, namely that in Eq. (36) 
we assumed the form of U of Eq. (25), thus assuming all the singular values of the matrix Q 
(22) to be nonzero, i.e. rank(S) = 4. In the rank(S) < 4-cascs it is still possible to write 
a decomposition of X of the form (46), but in these cases we have even more freedom in 
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choosing P\FP 2 (and also in choosing ft, though this freedom is not important). We now 
consider the different ranks of X seperately. 

4.2.0. rank(S) = 

The case rank (£) = is trivial, since in this case C» = and hence also C b = & = 0. 

4.2.1. rank(S) = 1 

For states where rank (£) = 1 we write S = diag (a, 0, 0, 0), so we have C" = a. From Eq. (22) 
we have Q = U*T,U^ and hence Qij = aU^U^. According to Eq. (18) a measurement in the 
standard basis will yield an average concurrence of 

k k 

so we have C b — C" for these states. We notice that in fact every measurement will do for 
these states. Since the standard basis works for all states, all measurement bases will work 
for any particular state. 

4.2.2. rank(S) = 2 

For the rank(S) = 2 case we apply the previously mentioned result by Walgate et al. [4]. 
Writing £ = diag (<ri, <J2, 0, 0) = U T QU as before we now have C" = <7i + oi- As discussed in 
Sec. , the column vectors of the unitary matrix U* give an optimal joint measurement basis. 
According to Walgate et al. we can write the first two (orthogonal) column vectors in the 
form 

u *i = c i\0)c\vo)d + c 2 |1)c|»?i)d , 

U* 2 - C3|0) c |^>d + c 4 |l)chf% , (49) 

with (r?o|?7o") = (v^Vi) = 0- A local measurement by C in the standard basis followed by a 
local measurement by D in either the {\rjo), \Vo)} or {l^i)) \Vi)} basis is equivalent to a 
joint measurement in the basis 

{\0)c\vo)d, \1)c\vi)d, |0> c |^)d, \1)c\Vi)d} ■ (50) 

The unitary matrix V, defined such that the column vectors of V* are these four vectors, is 
indeed of the form (26), and hence the average concurrence after this measurement is 

C = ^\(V T QV)kk\=E\( W ^ W ')kk\ . 
k k 

where we have defined the matrix W = V^U. From Eq. (49) we see that the first two column 
vectors of W* are W* x = (ci,c 2 ,0, 0) T and W* 2 = (0, 0, c 3 , c 4 ) T , so the average concurrence 
after this measurement is 

C = N^fei) 2 + <r 2 {W* k2 f\ = (j\ (| Cl | 2 + |c 2 | 2 ) + a 2 (|c 3 | 2 + |c 4 | 2 ) =*i+<T2, (52) 
k 
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attaining the upper bound C", and hence C b = C" also for these states. Even though classical 
communication between the assistant parties is needed in this scheme, we know from Sec. 
that it is possible to find another local measurement basis where classical communication is 
not necessary. 

4.2.3. rank(S) = 3 

For rank(S) = 3 the extra freedom in the choice of U of Eq. (36) is only 



U — >U 



( 1 \ 
10 
10 

\o o o J 



(53) 



but since we require det(£7) = 1, we have no more freedom in this case than in the general 
case, and hence the results for the rank(S) = 4 states apply here as well. So here we have 
C b < C» in general. 

5. Numerical Results 

In the previous section we have identified the set of states for which the concurrence of 
assistance can be achieved by local von Neumann measurements. We found that even though 
there was quite a lot of freedom in choosing these states, the set is of measure zero on the full 
state space, so for a randomly chosen state C" > C b . 

It would be interesting to see more quantitatively how big the advantage of joint measure- 
ments is. So far we have seen that there exist examples for all the extreme cases C b — C = 1, 
C b = C" = 0, and C b = 0,C = 1. But how large can we expect the advantage of joint 
measurements to be for an arbitrary state? 

To answer this question we have calculated the two quantities, C" and C , for a large 
number (~ 10 6 ) of arbitrary states. Our earlier analysis provides us with a numerical method 
to find C b for a given state, Eq. (35), and the concurrence of assistance is easily found from the 
formula (12). Picking pure states of four qubits arbitrarily is done by choosing the individual 
components Eq. (13), as independent random gaussian variables with mean value zero, 
each multiplied with a uniformly distributed phase factor [16]. The joint distribution function 
of these is indeed spherically symmetric, and therefore, when we have normalized the states, 
they will give a uniform distribution on the unit sphere in C 16 , which is what we are looking 
for. 

The results are presented in Fig. 2, showing the distribution of 25000 arbitrarily chosen 
states plotted according to their values of C" and C b (plotting all 10 6 states doesn't change 
the appearance of this figure). Here we see a clear accumulation of points close to the limiting 
C" = C b line, indicating that the gain of using joint measurements is in fact typically quite 
small. This is further substantiated by figure 3, showing the distribution function P(x) for 
the relative gain, x = (C ( - C b )/C b , normalized so that / P(x)dx = 1. Here we see that P(x) 
is decreasing rapidly (roughly exponentially) for larger x and the average gain is calculated 
to be only 

— )=4.6%. (54) 
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Fig. 2. The distribution of states according to C" and & , (25000 trial states). 




Fig. 3. The distribution function for the relative gain, (C" - C b )/C b , (10 6 trial states). 
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Fig. 4. The distribution functions for & (solid) and C b (dashed), (10 6 trial states). 

For a comparison of the two functions and C b , we have also plotted their (normalized) 
distribution functions separately, Fig. 4. The mean values of these distributions are found to 
be 

(C J ) = 0.778 and (C b ) = 0.745 . (55) 




0.2 - 



I ' ' ' ' 1 

0.2 0.4 0.6 0.8 1 

Joint Concurrence 

Fig. 5. The distribution of states according to the average of C* and C b over the six possible pairs 
of keepers of the final state, (25000 trial states). 

So far the problem we have considered has been for the assistant parties, C and D, to make 
A and B end up with as much entanglement (concurrence) as possible. Another interesting 
quantity would be the average of C" and C b over the six possible pairs of keepers of the 
final state. The figures 5, 6 show the distributions of these quantities. The mean values of 
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Fig. 6. The distribution functions for the average of C" (solid) and C (dashed) over the six possible 
pairs of keepers of the final state, (10 6 trial states). 



the two distributions in Fig. 6 will necessarily be the same as those in Eq. (55) since the 
states are chosen at random, but we notice that the variances are considerably smaller. This 
indicates that states with high C" or C b for one pair tend to have low C" or C b for some 
of the other pairs and vice versa. We stress that these are the results for arbitrary states. 
However, we can find examples of states violating the general tendencies. One would be the 
state \ip) = |0000) + |HH), for which local measurements by any two parties can produce a 
maximally entangled state shared by the other two, and another would be a separable state 
like \ip) = 1 0000), for which it is not possible for any parties to produce any entanglement at 
all. 

5.1. Nonprojective measurements 

An interesting extension of the current work would be to investigate the effect of nonprojective 
POVMs. If we allow joint measurements we have seen that von Neumann measurements are 
sufficient, see the discussion following Eq. (24), but are these also sufficient to achieve the 
concurrence of local assistance, C^l 

At most one of the assistant parties (the one who measures first) needs to do a nonpro- 
jective measurement. For after the first measurement the total state is a pure 3-party state 
and the problem is reduced to that of finding the concurrence of assistance for this state. For 
this we know from Sec. that a von Neumann measurement is sufficient. 

We have searched among states and measurement schemes where the first measuring party 
(C) applies a four-outcome POVM, and we have indeed found examples of states where a 
nonprojective measurement would yield more average concurrence than the best von Neumann 
measurement on the same state, i.e. > C b . A particularly simple example is the state 

|V>)abcd = -^(| $+ )ab| + +)cd + |3»~}ab|00) C d + |*~)ab|11)cd) , (56) 
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for which & = 0.8801 whereas & > 0.8978, i.e. an increase of about 2.0% (for this state 
C" = 0.9205). The inequality sign stems from the fact that we have only been considering four- 
outcome POVMs. The largest difference wc have found is for a randomly generated state, for 
which we saw an increase of 3.5%. For random states the increase is typically several orders 
of magnitude smaller than this. For certain states we have found that the maximum average 
concurrence is different if D measures first, indicating that classical communication between 
the two assistant parties might be necessary in doing the optimal nonprojective measurement. 

6. Discussion and Open Questions 

We have studied the concurrence of assistance for 4-qubit states where 2 parties are assist- 
ing. We introduced three versions of the concurrence of assistance, namely & where the 
assistants are allowed to perform joint measurements, & where they are restricted to local 
measurements, and CP where they are restricted to local von Neumann measurements. 

For joint measurements we have obtained the simple formula Eq. (24) for the concurrence 
of assistance, C". For local von Neumann measurements we have found that classical com- 
munication of measurement outcomes between the two assistant parties is not necessary, but 
for four-outcome POVMs we have seen that communication does help. A main result is the 
decomposition of Eq. (46), telling us (i) the concurrence of assistance, (ii) whether there is a 
local von Neumann measurement that can be used to achieve this and (iii) what that local 
measurement would be. By using numerical methods we have seen that even though we can 
construct states for which the gain of concurrence in using joint measurements is anywhere 
between and 1, for the typical state C" is only around 5% higher than C b . Finally we have 
seen that in general nonprojective measurements can do a little better than von Neumann 
measurements, i.e. the inequality > & is strict. 

Several open questions arise naturally from the present work. The criterion of Sec. 
to determine the set of states for which local measurements are sufficient to achieve the 
concurrence of assistance only takes into account von Neumann measurements. Whether 
there are states for which £P < = C" , and whether there are states for which C" > = C b 
are both open questions. Other interesting problems would be to see whether communication 
between the assistant parties is needed for some of the states where C b < C^. Finally, an 
obvious question is what the results are if we use the entanglement (5) as our measure instead 
of the concurrence. 
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Appendix A: Proof of Eq. (23) 

This is a generalization of a result in Ref. [14] (p. 430). 

Let Q be a square matrix with singular values <Jk , and write its Takagi decomposition (22) 
as U T QU = S. Let V represent an arbitrary right unitary matrix, VV^ = 1. Then 



max |Tr (V T QV) | = max |Tr ((C/ t U) T S(C/ t U)) | 

= max |Tr (V T EU) | = max |Tr (Y,VV t ) 



= max 
v 



°k (VV T ) 



kk 



<max]>>|(VU T ) fefc |<5>, (A.l) 



where the last inequality follows from the fact that 



VV 1 



'kk I 



<Ei^i 2 = ( yt/t )^ = 1 - 



Now assume that there exist a right unitary matrix Vb satisfying 

EI(^ T QVoU>]>>- 



k k 



(A.2) 



(A.3) 



Let its diagonal entries be (V QVo)kk — r^t ' k with real and positive, and define the 
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matrix V\ = VqD, where D is diagonal and Dkk = e~* 5fe / 2 . Since V\ is right unitary, and 

|Tr (V^QV 1 ) |=5>>E**' (A.4) 

k k 

in contradiction with Eq. (A.l), we conclude that 

EK^OJ^E^ ( A - 5 ) 

k k 

for all right unitary matrices V. 
Appendix B: Proof of Eq. (34) 

Let U be unitary. Write the real and imaginary parts of U and find the singular value 
decompositions of the real part, 

u = u R + iu I = 6 1 v R dl + iUj, (b.i) 

where 0\ and 2 are both real and orthogonal. Now look at the matrix U — 0\UO^ = 
Y>r + iA. Unitarity, UU* = &U=1 implies 

E| + AA T = Y, 2 R + A T A = 1 and AY, R - H R A T = , (B.2) 

so the real matrix A is normal, AA T = A T A, and hence it has an eigenvalue decomposition 
[14], 

A = VAV* , (B.3) 

with V unitary. The matrix 

AA T = AA^ = KAVtVA'Vt = V |A| 2 (B.4) 

is diagonal according to Eq. (B.2). If all its entries are distinct, V can only be making 
permutations of these, and hence it can also only do so in Eq. (B.3), meaning that A is 
diagonal. 

If some of the diagonal entries of AA T are identical, we first apply the real and orthogonal 
operator P permuting the diagonal entries, to ensure that the identical ones appear last on 
the diagonal of AA T . Thus we define U = PUP T = S fi + iA. 

For the special case of two identical and two distinct entries, the matrices are now 

iI T = diag(|A 1 | 2 ,|A 2 | 2 ,|A 3 | 2 ,|A 3 | 2 ) , 

± 2 R - diag (l - |A!| 2 , 1 - |A 2 | 2 , 1 - |A 3 | 2 , 1 - |A 3 | 2 ) , 

A = diag (Ai, A 2 ) ® A' . (B.5) 

From the unitarity condition (B.2), AY. R — T^ R A T or Aij&j — &iAji, we see that the submatrix 
A' of Eq. (B.5) is symmetric, so we can write its eigenvalue decomposition A' = V'A'V' T with 
V real and orthogonal and with A' = diag (A 3 , A4). Defining the matrix V = 10 V we finally 
get 

VPO x UOlP T V T = vllV T = D , (B.6) 
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where D = S fl + idiag (Ai, A 2 , A' 3 , A4). 

This generalizes straightforwardly to other degeneracies, and hence we have the desired 
decomposition (34), where we identify Oi = VPO\ and O2 = VPOi- 

This decomposition is not completely unique. The eigenvalues and singular values that 
we have found are unique, so the matrix D is unique up to certain permutations between 
the real values and between the imaginary values of its entries. For D to remain unitary 
all the diagonal entries must still be of the form e ltp , and hence the only allowed form of 
permutations can be Im (Da) <-> Im (Djj) with Im (Djj) = — Im (Da). Since Re (Da) > per 
construction, this is the same as interchanging Da <-> Djj, and hence we conclude that D is 
uniquely defined up to permutations of the diagonal entries, which we shall refer to as the 
phases of U, e lSi , defined such that — ^ < 5 j < §. 

For a given ordering of the phases in D, the 0\ and O2 are unique up to simultaneous 
transformations of the kind 0\ — > 0\0, 2 — > 2 0, where O is a real, orthogonal matrix 
leaving D invariant, ODO T = D. In the case where all the phases are different, O is a 
diagonal matrix with ±1 on the diagonal. Otherwise it is a block diagonal matrix, with the 
sizes of the blocks determined by the multiplicity of the phases. 



